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ABSTRACT 

The Baryon Acoustic Oscillations (BAO) are features in the matter power spectrum on scales of 
order 100 - 150 /i~^Mpc that promise to be a powerful tool to constrain and test cosmological 
models. The BAO have attracted such attention that future upcoming surveys have been designed 
with the BAO at the forefront of the primary science goals. Recent studies have advocated the 
use of a spherical-Fourier Bessel (SFB) expansion for future wide field surveys that cover both 
wide and deep regions of the sky necessitating the simultaneous treatment of the spherical sky 
geometry as well as the extended radial coverage. In this paper we present and extended analysis 
of the BAO's using the SFB formalism by taking into account the role of non-linearities in the 
oscillations observed in the galaxy power spectrum. The SFB power spectrum has both radial and 
tangential dependence and it has been shown that in the limit that we approach a deep survey the 
SFB power spectrum is purely radial and collapses to the Cartesian Fourier power spectrum. This 
radialisation of information is shown to hold even in the presence of redshift space distortions 
and 1-loop corrections to the galaxy power spectrum albeit with modified tangential and radial 
dependence. As per previous studies we find that the introduction of non-linearities leads to a 
damping of the oscillations in the matter power spectrum. 

Key words: : Cosmology- Cosmic Microwave Background Radiation- Large-Scale Structure of 
Universe - Methods: analytical, statistical, numerical 



1 INTRODUCTION 

Observations of the cosmic microwave background (CMB) and large- 
scale structure (LSS) will carry complementary cosmological infor- 
mation. While all-sky CMB observations, such as NASA'a WMAP or 
ESA'a Planck experiments, primarily probe the distribution of matter 
and radiation at redshift z = 1300, large scale surveys such as ESA's 
Euclid or the SKA will provide a window at lower redshifts on order 
z ^ — 2. The study of large scale structure appears to be a promis- 
ing candidate in the study of the influence and role of the dark sec- 
tors in the standard model of cosmology. One particular phenomena 
of interest are the Baryon Acoustic Oscillations (BAOs) that manifest 
themselves in the matter power spectrum of galaxy clusters on cosmo- 
logical scales of order 100h~^ Mpc. These oscillations in the matter 
power spectrum are generated just before recombination through the 
interplay between a coupled photon-baryon fluid and gravitationally 
intera cting dark matter ([Sunyaev and Zeldovich 1970, Peebles and Yu| 
[T970| [Eisenstein et al|2005[|Seo and Eisensteinl2003]|2007l L 



The scale of the peaks and oscillatory features of the BAOs 
promises to be an important cosmological tool that acts as a standard 
ruler from which we can investigate and constrain dark energy pa- 
rameters (see Eisenstein et al[p005] ); [Amendola, Quercellini and Gi 



'allongo ('2005]); [Dolney, Jain and T^ada| ( |2006l ); [Wang and Mukher 
jee ( 2006 ) for a small selection or representative literature), neutrino 
masses (Goobar, Hannestad, Mrtsell and Tu 20()6|, modified theories 



of gravitation ( ^Alam and Sahni,2006^^Lazkoz, Maartens and Majerotto 
2006[) and deviations from the stand ard model of cosmology ([Garcia 



Bellido and Haugboell e[2008l [20091 [February, Clarkson and Maartens 
2012[ )). Significant attention has been devoted to the BAOs and they 



were first detected with SDSS data ( [Eisenstein et al| 2005", AdelmanJ 
McCarthy et al|2008] ) and in subsequent surveys ( [Colless et al,2003^ 



Percival et al|2007l ). 



The BAOs have been studied using standard Fourier space de- 
Eisenstein[[20 0? "2007), real space analysis 



compositions (|Seo anc 
( [Eisenstein et al[[2005'^ 



[^ [20T0l 



Juszkiewicz, 



[Slosar, Ho, W hite and Louis 20091 [Xu et[ 



iellwing and van de Weygaert||2012| ), in 2D 
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spherical harmonics defined on thin spherical shells ([Dolney, Jain and 



Takada 2006 j, but al so in a spherical Fourier-Bessel expansion ( [Rassat 



and Refregier'20121). It is important to note that different frameworks 



will make use of different information and will therefore have dif- 
ferent constraining power for different cosmological parameters em- 
phasising the complementarity of mixed studies (Rassat et al 2008 ). 
Previous studies, having predominantly focused on projected 2D sur- 
veys, have discarded statistical information by projecting galaxy posi- 
tions into tomographic redshift bins however, such a loss of informa- 
tion could be avoided by adopting a full 3D description ( [Asorey et al| 
[2012| ). 

Upcoming large scale structure surveys will provide cover for 
both large and deep areas of the sky and this will necessitate a for- 
malism that can provide a simultaneous treatment of both the spher- 
ical sky geometry as well as an extended radial coverage. A natural 
basis for such a survey is provided by the spherical Fourier-Bessel 
(SFB) decomposition, see ([Heavens and Taylor||1995 [Fisher et"aT 



1995 , Percival et al TOW, Cast ro, Heavens an d Kitchmg"2005; "Er 
dogdu et al 2006, Abramo, Reimberg and Xavier 2010, Leistedt et 
al|2011 Shapiro, Crittenden and Percival 2011 , Rassat and Refregier 



2012[|Lanusse, Rassat and Starck,20 12^, Asorey et al,2012 ) for an in- 



complete selection of literature on the subject. In this prescription we 
expand a 3D tracer field, such as the galaxy density contrast, using 
the radial {k) and tangential (i.e. along the surface of a sphere) (/) 
dependence. 

The galaxy matter power spectrum, as an initial approach, is of- 
ten modelled using cosmological perturbation theory (PT) at linear or- 
der. However, linear PT will break down due to the nonlinear growth 
of structure under gravitational instability. This necessitates the use of 
non-linear (or higher order) PT. In addition real world observations 
are plagued by the presence of nonlinear contributions to the power 
spectrum. The main sources of non-linearity that we must take into 
account include: 

• Non-linear evolution of matter density field and it's perturba- 
tions. This will give rise to correction to the linear matter power spec- 
trum. 

• Non-linear galaxy clustering bias arising from a non-linear map- 
ping between the underlying matter density field and observed col- 
lapsed objects (e.g. galaxies or DM Haloes). Galaxy clustering bias 
also modifies the matter power spectrum. 

• Redshift space distortions arising from the peculiar velocities of 
galaxies manifest themselves in the observed redshifts of galaxy sur- 
veys. RSD give rise to systematic effects in the redshift space matter 
power spectrum that differ from the real space matter power spectrum. 

• Non-Gaussian corrections to the initial density field. In this paper 
we assume Gaussian initial conditions. 



In this paper we will follow the construction outlined in [Rassat] 
[and Refregier ( 2012 ) and generalise the method to study the role of 
redshift space distortions (RSD) and the non-linear (NL) evolution of 
density perturbations. Previous investigations have used standard per- 
turbation theory (SPT), galaxy bias models and Lagrangian perturba- 
tion theory to characterise the role of various non-linear corrections to 
the BAO signal using the 3D Fourier power spectrum P{k) ( [Jeong and| 



Komatsu |2006i|Nishimichi et al|200 7[[Jeong and Kom atsu|2009l|Na 
mura, Y amamoto and Nishimichi 20081 [Nomura, Ya mamoto, Huetsi 
and Nishimichi,2009j . These nonlinear corrections can be reassessed 



within the SFB framework to aid our understanding of how real world 
effects can impact the radialisation of information. 

Recent work ( |Asorey et al 2012 ) utilising the SFB formalism has 
focused on how to recover the full 3D clustering information including 
RSD from 2D tomography using the angular auto and cross spectra of 
different redshift bins. Traditionally, RSD measurements have been 
made through spectroscopic redshift surveys such as the 2dF Galaxy 
Redshift Survey ( ^Colless et al 2003 ) and the Sloan Digital Sky Survey 
( [York et al|2000[ ) with photometric surveys often being neglected be- 
cause of the loss of RSD through photometric redshift errors. Upcom- 
ing surveys, spectroscopic and photometric, such as the Dark Energy 
Survey (DES www.darkenergysurvey.org), Euclid (sci.esa.int/euclid). 
Square Kilometer Array (SKA www.skatelescope.org). Physics of 
the Accelerating Universe Survey (PAU www.pausurvey.org) ( [Bentez[ 
[et al|[2009| ). Large Synoptic Survey Telescope (LSST www.lsst.org) 
or the Panoramic Survey Telescope and Rapid Response System 
(PanStarrs pan-starrs.ifa.hawaii.edu) offer the possibility of investi- 
gating the BAO and RSD through angular or projected clustering 
measurements Bentez et"aT| (2009); Nock, Perc ival and Ross ( |2010| ); 
[Crocce, Fosalba, Casta nder and Gaztanaga ( 2d T0|); Gaztanaga et al 



([201 1[); Laureij s et al [ ( ^20 1 1 ) ; Ross, Percival, Crocce, Cabr and Gaz- 
[tanaga| ( 20TT] L 



As RSD and distortions arising from an incorrect assumption for 
the underlying geometry are similar ( [Alcock and Paczynski|1979] ) the 
analyses of RSD using 3D data has to be used in conjunction with ge- 
ometrical constraints. Approaches based purely on angular correlation 
functions, depending only on the redshift and angles, do not depend on 
the background cosmological model making angular clustering mea- 
sures considerably simpler in this respect. The SFB is something of a 
mid-point between these two extremes and will, in general, be sensi- 
tive to the choice of fiducial concordance cosmology. 

This paper is organised as follows. In ^we discuss the Spher- 
ical Fourier-Bessel (SFB) Expansion in the context of redshift- space 
distortions. In ^we add the perturbative corrections to linear power- 
spectrum. Results are discussed in ^and finally the conclusions are 
presented in ^ Discussions about finite size of the survey and dis- 
crete Spherical Fourier-Bessel (SFB) transforms are relegated to an 
appendix. 

Throughout we will ad opt the WMAP 7 cosmological param- 



eters |Komatsu et al 



2011 



h -- 

0.112, Qa = 0.725, erg = 0.816. 



0.7, Qth' 



0.0226, Qc/i^ 



2 REDSHIFT SPACE DISTORTIONS 

2.1 Spherical Fourier-Bessel (SFB) Expansion 

Spherical coordinates are a natural choice for the analysis of cosmo- 
logical data as they can, by an appropriate choice of basis, be used to 
place an observer at the origin of the analysis. Upcoming wide-field 
BAO surveys will provide both large and deep coverage of the sky and 
we therefore require a simultaneous treatment of the extended radial 
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coverage and spherical sky geometry. For this problem, the spherical 
Fourier-Bessel expansion is a natural basis for the analysis of random 
fields in such a survey. 

We introduce a homogeneous 3D random field ^(H, r) where Q 
defines the position on the surface of a sphere and r denotes the co- 
moving radial distance. The eigenfunctions of the Laplacian operators 
are constructed from products of the spherical Bessel functions of the 
first kind ji (kr) and spherical harmonics Yim {^) with eigenvalues of 
—k^ for a 2-sphere. Assuming a flat geometry, the random field can 
be harmonically decomposed into the spherical Fourier-Bessel basis, 
which is both complete and orthonormal, with the expansion given by: 




and the corresponding inverse relation given by: 

^imik) = y d\^{r)kji{kr)YiU^). (2) 

In our notation {Im} are quantum numbers and k is the wavenum- 
ber. Note that the 3D harmonic coefficients, "^imik) are a function of 
the radial wavenumber k. This decomposition can be viewed as the 
spherical polar analogy to the conventional Cartesian Fourier decom- 
position defined by: 

*W=(^/rf'fc*We"^ (3) 
*(k)=(^/rf'-*(r)e-*-. (4) 

The Fourier power spectrum, P^^, is defined as the 2-point correla- 
tion function of the Fourier coefficients ^(/c), 

(^(k)^*(kO) = {27vfP^^{k)6' (k - k^) . (5) 

Similarly we can define a 3D SFB power spectrum, Ci{k), of our 
random field by calculating the 2-point correlation function of the 3D 
harmonic coefficients, 

(*,„(fc)vl/;,„,(fc')) = C,(fc)5^°(fc - fc')4"5w- (6) 

It is possible to relate the Fourier coefficients ^(k) with their SFB 
analog ^iyn{k) through the following expression, 

^^Imik) = J^^^^ j dVLk^{\<i)Ylm{^k) (7) 

where the angular position of the wave vector k in Fourier space is 
denoted by the unit vector tt{Oi^(j)k). The Rayleigh-expansion of a 
plane wave is particularly useful in connecting the spherical harmonic 
description with the 3D Cartesian expression. The second expression 
we present here is derived by differentiating the first and will be used 
in the derivation of redshift space distortions: 

e^k-^ = 47r^i^ji(kr)Yin,(nk)Yin,(Q); (8) 

Im 

• n)e'''-" = 47r^i^j;(kr)Yi„,(Qk)Yi„,(Q). (9) 

Im 



In general the radial eigenfunctions are ultra- spherical Bessel func- 
tions but they can be approximated by spherical Bessel functions when 
the curvature, of the Universe, is small. Throughout this paper we will 
use j'l (x) and ][' (x) to denote the first and second derivatives of ji (x) 
with respect to its argument x. The expressions for the first and second 
derivatives are given in Eq.(B2) and Eq.(B3). Imposing a finite bound- 
ary condition on the radial direction will result in a discreet sampling 
of the k-modes. This will be discussed in more detail later. 

2.2 Selection Functions and Window Functions 

In order to consider realistic cosmological random fields, such as the 
galaxy density contrast, we need to take into account the partial obser- 
vation effects arising from finite survey volumes. This discussion will 
follow that of jRassat and Refregier| ( |20l"2| ), including the discussion 
of the radial and tangential limits of the selection function. The selec- 
tion function denotes the probability of including a galaxy within the 
survey. 

An observed random field (r) can be related to our underly- 
ing 3D random field by a survey-dependent radial selection function 

^°^^(r) = 0(r)^(r). (10) 

It is possible to introduce and analogous tangential selection function 
but we will, as per jRassat and Refregier| ( |2012] ), neglect this possibility 
assuming that we have data for the full sky. The selection function is 
typically normalised by the characteristic survey volume, V, defined 
such that as we take the volume to infinity we recover the underlying 
random field, 

lim ^"^'(r) = ^(r). (11) 

As homogeneity is broken, due to the introduction of the selection 
function, the SFB power spectrum now depends on /c, k' and / and is 
given by: 

(fc)*?'t'> = Cf^k, k')5u.5^^,, (12) 

but can be shown to be expressed as follows (e.g. |Rassat and Refregier| 
( [20T2I )): 

C?0(A;i,A:2)= 0^ j k'^dk' I^{ki,k')I^{k2,k')Ps5{k') (13) 

where the window function is defined as: 

l[^\k,k')^ j drr'^(t){r)kji{kr)ji{k'r). (14) 

The SFB power spectrum tends to rapidly decay as we move away 
from the diagonal k — k' and it will often be much more useful to 
focus purely on the diagonal contribution. 

2.2.7 Gaussian Selection Function 

A commonly adopted selection function is the Gaussian function, 

0(r) =e-(^/^°^', (15) 
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such that ro is a radius parameter defined by th e normahsation of the 
survey volume, V — vr^/^ro. It can be shown ( Rassat and Refregier 
|2012| l that the window function may be integrated analytically and 
reduces to: 



y exp 



-^0 



rlkk' 



(16) 



2.2.2 Radial Limit 

The radial limit corresponds to full radial coverage such that, 

0(r) ^ 1 Vr (17) 

and the Gaussian window function is significantly simplified to 

Ii{k,k') = ^5{k-k'). (18) 

Using this asymptotic window function the power spectrum reduces 
to the Cartesian Fourier power spectrum ( |Rassat and Refregier|2012| ): 



Cf'^'ik^k') ^ Ci(k)5{k - k') ^ P{k)5{k - k'). 



(19) 



2.2.3 Tangential Limit 

The second case is defined in the limit that the selection function 
probes only a thin shell of the field at a distance r* ( ,Rassat and Re-| 
|fregier|2012| ). 



(j)(r) — r^S (r — r*) . 



(20) 



normalised such that V = 47rrJ. The corresponding convolution ker- 
nel, neglecting redshift space distortions, reduces to. 



ll''\k,k')=rlji{kr.)ji{k'r.) 



2.3 Power Spectrum Evolution and Bias 



(21) 



The linear matter power spectrum at an arbitrary redshift z can be 
calculated by evolving the power spectrum according to 



Pssik.z) = D\z)I\ni{k)T\k) 



(22) 



where T(k) is the conventional transfer function. The linear growth 
factor D{z) is normalised to unity at present times and allows us to 
quantify the change in the linear growth of matter with respect to red- 
shift. In a realistic galaxy survey, observations of galaxy density field 
contrast 6g act as a tracer of the underlying matter distribution 6m ' 



Sg = b{z, k)Sn 



(23) 



such that b{k, z) is the galaxy bias. In general, the bias may be both 
redshift and scale dependent but we do not investigate the role of bias 
to any significant detail in this paper instead adopting a common ad- 
hoc functional form for the bias, b{z) = \/l + z. Throughout we 



will assume that the bias is scale independence and therefore only 
modulates the overall amplitude of the galaxy power spectrum Pg (k) 



)D\z)P,r.i{k)T\k). 



(24) 



2.4 Redshift Space Distortions 



The measured distribution of galaxies is not without limits though as 
various systematic and survey dependent errors become more impor- 
tant. In practice, the observed galaxy redshift distributions are dis- 
torted due to the peculiar velocity of each galaxy. The anisotropy 
generated by the peculiar velocities is known as redshift space dis- 
tortion. Although this distortion of the measured redshifts will neces- 
sarily complicate the cosmological interpretation of the spectroscopic 
galaxy surveys, the redshift space distortions are currently one of the 
most optimistic probes for the measurement of the growth rate of 
structure formation and, as a result, an interesting probe of models 
for dark energy and modified theories of gravity. 

The effect of the redshift space distortions on the matter power 
spectrum can be split into two effects, the Kaiser effect and the Finger- 
of-God effect. The Kaiser effect corresponds to the coherent distortion 
of the peculiar velocity along the line of sight with an amplitude con- 
trolled by the growth-rate parameter. The Finger-of-God effect arises 
due to the random distribution of peculiar velocities leading to an in- 
coherent contribution in which dephasing occurs and the clustering 
amplitude is suppressed. It is thought that the suppression of the am- 
plitude is particularly important around the size of halo forming re- 
gions. 

For an isotropic structure in linear theory, the Kaiser effect means 
that an observer will measure more power in the radial direction than 
in the transverse modes. The amplitude of this distortion is modulated 
by the distortion parameter 



/3 = 



1 d\nD{a) 



(25) 



b{z) b{z) d\na b{z) 

meaning that the RSD can be used as a probe of the growth of structure 
and of the galaxy clustering bias function b{z). 

The effect a peculiar velocity, or a departure from the Hubble 
flow, v{r) at r is to introduce a distortion to the galaxy positions in 
the redshift space s: 

s(r) = r + v(r) • Q. (26) 

We denote the harmonics of a field ^(r) when convolved with a se- 
lection function, (/)(s), by ^im{k). These harmonics take into account 
the redshift space distortions: 



-^Imik) 



(iQ0(s)^(r)j,(A:s)y;^(Q). (27) 



The Fourier transform of the linearised Euler equation can be used to 
relate the Fourier transform of the density contrast, 5 (k), to that of the 
peculiar velocity field v{y): 

/(k) 



v(k) = -i/3k- 



A;2 



(28) 
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where b is the Unear bias parameter. Following the procedure outline 
in [Heavens and Taylor] ( [1995] ), we can establish a series expansion 
in P such that the lowest order coefficients ^[^(k) are obtained by 
neglecting the redshift space distortions: 



,f,(l)/ 



/ o POO 

Mnl(k) = J^ k'dk'^ 
V TT Jo 



2 

^ Jo 



l['\k\k). 



(29) 
(30) 

(31) 



The kernels {k\ k) and {k\ k) define the convolution and are 
dependent on the choice of selection function. Note that Ii^\k\ k) 
is simply the window function we encountered previously in Eq.(T4|. 
The kernels can be shown to be: 



r(0) 



ll'\Kk') 



{k^k') — j dr (j){r) k ji{kr) ji{k' r) 



(32) 



drr'^ k ^ {(t){r)ji{kr)) j[{k'r). (33) 



The lowest order corrections due to redshift space distortions are 
therefore encapsulated in {k). We can define a set of power spec- 
tra by using these harmonic coefficients: 



(*r™(fc)*f4,(fc')> 

(*L(fc)*f4'(fc')> 



C\""\k,k')5^D{k-k')5u'5„ 



Cr'{k,k')6u, 



(34) 
(35) 



We can construct a generalised power spectrum by using the common 
structure between Eq.pO^ and Eq.jSTJ: 

C\''^\k,M)^ (^)V ^'^dk' l[''\k,,k')l\^\k2.k')Pss{k'). 

(36) 

The total redshifted power spectrum will be given by a sum of the 
various contributions: 

Cl (ki , /C2) = (/Cl , /C2) + 2 (ki , /C2) + C[''^ (ki , /C2) . (37) 

If we ignore the effects introduced by the selection function, i.e. set 



(/)(r) = 1, then we recover the result of Castro, Heavens and Kitching 



( |2005| ) for the unredshifted contributions: 

C'i°°'>ik,k) = Pss{k). 



(38) 



These expressions hold for surveys with all-sky coverage. In the pres- 
ence of homogeneity and isotropy the 3D power spectrum will be in- 
dependent of radial wave number /. The introduction of a sky mask 
breaks isotropy and introduces additional mode-mode couplings, the 
analysis will be generalised to this case in the next section. In the 
above equations we neglect a number of additional non-linear terms 
including General Relativistic corrections, velocity terms and lensing 
terms. It is also possible to adopt a full non-linear approach to RSD 
where the non-linear spectrum has significantly more complicated an- 
gular structure than in linear theory [Shaw and Lewis| ( |2008| . The RSD 
information will be dependent on the relative clustering amplitude of 
the transverse modes and the radial modes [Asorey et allp012| ). Our 



ability to recover information and the extent to which the information 
radialises will naturally depend on the geometry of the survey and 
which modes we are able to include. 



2.5 Realistic Surveys 

The results that have been discussed above are somewhat idealised in 
the sense that we assume all-sky coverage with no noise. In realistic 
surveys we will often need to take into account the presence of a mask 
(relating to partial sky-coverage) and noise. If the noise is inhomoge- 
neous we will be presented with a further complication. For partial sky 
coverage we find mode-mode couplings in the harmonic domain that 
result in the individual masked harmonics being described by a lin- 
ear combination of our idealised all-sky harmonics. We do not discuss 
the role of partial sky-coverage in much detail but do present results 
generalising our formalism to include a survey mask. 



2.5.7 Partial-Sky Coverage and Mode Mixing 

Large scale surveys do not, generally, have full- sky coverage. Instead 
the information regarding sky-coverage is encapsulated in a mask 
x(ft) which is unity for areas covered in the survey and zero for re- 
gions outside the survey. The field harmonics are therefore modulated 
in the presence of a mask: 

*i™(fc) = j s^ds JdCl [<k{s)x{n)] ^{r)ji{ks)YtU^). 

(39) 

The convolved power-spectra in the presence of the mask takes the 
following form: 

c[''^\ki,k2) = (^)'ee/ ^'^^'/ 

X Wtj{k,,k')w[fj{k2,k'')j^^^ (40) 



Cf {XlmXlm) 



(41) 



where I11I2I3 is the Gaunt integral (see Eq.( |B4| )). The convolved 
power spectrum is a linear combination of all- sky spectra and depends 
on the power spectrum of the adopted mask (see Appendix-|D]for de- 
tailed derivations.). 



2.5.2 Error Estimate 

The signal to noise for individual modes for a given power- spectrum 
can be expressed as: 



SCijk^k) 
Ci{k,k) 



21 + 1 



1 + 



1 



nCi{k, k) 



(42) 



Where n is the average number density of galaxies and the second 
term represents the leading order shot-noise contribution. For our re- 
sults we take n = 10"^/i^Mpc~^. 
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2.6 BAO Wiggles Only 

The BAOs can be isolated by constructing a ratio between the ob- 
served matter power spectrum Pfsi^) ^ theoretical matter power 
spectrum Pss{k) constructed from a zero-baryon transfer function in 
which the oscillations do not show up ( [Eisenstein et al|2005] ). Using 
these two power spectra, the ratio R^{k) will reduce the dynamic 
range and isolates the oscillatory features of the BAOs 



Wide and Shallow Survey: ro=100 /T^Mpc, / = 5 



(43) 



This ratio is clearly defined for the Fourier space power spectrum 
but an appropriate generalisation to the SFB formalism may be con- 
structed by calculating the ratio of the angular power spectra defined 



in Equation (13 1, with the matter power spectrum Ci (k) to the an- 
gular power spectrum with the zero-Baryon power spectrum CJ^^ (k) 
( |Rassat and Refregier|2012| ): 



RC.^.Clikl 



(44) 



It is important to note that the characterisation method (i.e. how 
we choose to construct our ratio) can affect the characteristic scale of 
the BAOs when we take into account non-linear effects. This means 
that care has to be taken when comparing results that implement dif- 
ferent methods ( [Rassat et al|2008| ). As an example we could construct 
our ratio by using the No- Wiggles transfer function of [Eisenstein et| 
[all ( |2005 ) or adopt an interpolation scheme to construct a smooth 
parametric curve [Blake et al| p006| ); [Percival et al] ( |2007| ); |Seo and| 
|Eisenstein|p007| ). A different choice of smoothed matter power spec- 
tra, cosmological parameters, growth history or similar can impact 
the phenomenological behaviour of the underlying physics (e.g. lo- 
cation of BAO peaks). Other methods for characterising the acoustic 
oscillation scales can be found, for example, in ( [Percival et al|2007{ 
[Nishimichi et al|2007| ). 



3 NON-LINEAR CORRECTIONS 

The role of nonlinear gravitational clustering can investigated in the 
spherical Fourier-Bessel formalism by incorporating higher order cor- 
rections to the power spectrum as described in perturbation theory. 
The approach we adopt here is standard perturbation theory (SPT) 
which provides a rigorous framework from which we can investigate 
the dependency of the radialisation of information on different param- 
eters in a fully analytic manner (IVishnia c[[1983| [Fry[[1984[ [Goroff, 
Grinstein, Rey and Wise 1986, Suto and S asaki[1991 [Makino, Sasaki 
and Suto 1992, Jain and Bertschinger 1994[ [Scoccimarro and Erie 
man, 1996 J . Standard perturbation theory is one of the most straight 
forward approaches to studies beyond linear theory and is based on 
a series solution to the hydrodynamical fluid equations in powers of 
an initial density or velocity field. The nonlinear clustering of matter 
arises from mode-mode couplings of density fluctuations and velocity 
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Figure 1. Slice in 1-space showing Rf (k) for = 5 (top-panels) and ^ = 50 
(bottom-panels) in a wide and shallow survey of ro = 100/i~^Mpc (left- 
panels) as well as for a deep survey of ro = 1400/i~^Mpc (right-panels). The 
solid-lines denote the C^^^^ term and the dashed-lines show the total power- 
spectra. 
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divergence as seen from the Fourier space equations. The role of per- 
turbation theory in the nonhnear evolution of the BAO in the power 
spectrum has been previously investigated (for an incomplete selec- 
tion of references please see:[fe)ng and Komatsu ( 2006 ); Nishimichi 
et al" ("2007"); "Nomura, Yamamoto and Nishimichi ( 2008 ); Nomura, 
Yamamoto, Huetsi and Nishimichi (2009 ); Taruya, Nishimichi, Saito 



and Hiramatsu| ( |2009| >; |Taruya, Nishimichi and Saito| ( |2010| )). In this 
paper we generalise these investigations to the SFB approach. 

Consider the hydrodynamic equations of motion for density per- 
turbations 5 such that our coming coordinates are denoted by x and 
the conformal time by 77: 

(5^(x, ry) + V • [(1 + ^(x, ry))v(x, ry)] = 0, (45) 

v'(x, T]) + [v(x, T]) ■ V] v(x, T]) + H(T])w(yi, T]) = -V0(x, ry), 

(46) 

V'(/)(x,7y) = ^7^'(^)^(x,ry), (47) 

where a prime denotes the derivative with respect to the conformal 
time and 1-L — a' /a. The rotational mode of the peculiar velocity v 
is a decaying solution in an expanding universe and can be neglected 
in this approach. We introduce a scalar field describing the velocity 
divergence: 



e(x,7y) = V • v(x,7y). 



(48) 



In our discussion we will focus on a description of the density pertur- 
bations and Fourier decompose the above equations to set up and solve 
a system of integro-differential equations. The Fourier decomposition 
of the perturbations are defined by: 



(49) 



e(x,ry) = I ^e(k,r;)e-^^- (50) 
The equations of motion can be decomposed as follows: 
(5'(x,r?) + e(k,r?) = (fki J d^fc2<5*^'(ki +k2 -k) (51) 



'^e(ki,,?)5(k2,r,), 

e'(k, n) + n{n)e{K n) + ln\r])S{K n) = 

,fc2(ki -ka 



(52) 



ZiKi-^ rh'2 



In order to solve these coupled integro-differential equations we intro- 
duce a perturbative expansion of our variables 

oo 

5(k,r,) = ^a"(r?)5„(k), (53) 

n=l 

oo 

e(k,,?) = ?i(,?)^a"(,?)G„(k). (54) 



The general n-th order solutions are given by: 

X Fn(qi, . . . ,qn)n7^i^i(qO, 



X C?n(qi,...,qn)nILi^i(q^), 



(55) 



(56) 



where the kernels Fn(qi, . . . , qn) and G'n(qi, . . . , qn) are given by 
I Jain and Bertschinger| ( |1994| ) 

, X v;^ Gm(qi, • • • ,qm) 

F4qi,...,qn)-2.^ (2n + 3)(n-l) 



k • ki 

(l + 2n)— ^ -i^n-m(qm + l, • • • , qn) 
^1 



k\ki •k2 



" Gn — m (qm + 1 5 



(57) 



C?n(qi,...,qn) ^ 



Cmlqi, 



(2n + 3)(n- 1) 

2 n — ra vQm + 1 5 • • • 5 Qn) 



+ n 



fc'(ki • ka] 

7.2 7.2 

^1^2 



Cn — m(qm + l, • • • , qn) 



(58) 



The kernel F„(qi, . . . , q„) is not symmetric under permutations of 
the argument qi . . . q„ and must be symmetrised 



^ Fn{m,...,Cln). 



(59) 



Permutations 

As an example the second order symmetrised solution is given by 

^2 



5 ^ 2 (ki • k2)^ ^ (ki • k2) / 1 



7 7 kfk 



2Z^2 
2 



+ 



+ 



1 



' 1.2 
^1 ^2 



(60) 



The corresponding second order matter power spectrum represents the 
linear matter power spectrum plus the additional higher order correc- 
tions. This calculation is made under the assumption that the first or- 
der density perturbations (k) constitute a Gaussian random field. 
The power spectrum up to second order is given by 



Pspt(A:, z) = D\z)Pi,n{k) + D\z)P2{k), 



(61) 



where Pun is the conventional linear matter power spectrum and the 
second order correction are given by 



P2(/C) =P22(/C) + 2P13(/C). 



(62) 



These terms correspond to the contributions to the 4-point correlation 
function from the (2,2)-order and the (l,3)-order cross-correlations. 
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The explicit form of these terms is given by (for full equations see ??) 

P22W = 2 j dVlin(|k-q|)[F|(q,k-q)]^ (63) 

Pi3(A;) = 3Piin(^) j d\Piin(^)F|(q,-q,k). (64) 

It should be noted that the analytical predictions arising from standard 
perturbation theory will eventually break down as the non-linear terms 
become dominant over the linear theory predictions. |Jeong and Ko-| 
[matsu ( 2006 ) demonstrated that one-loop standard perturbation theory 
was able to fit N-body simulations to greater than 1% accuracy when 
the maximum wave number ki% satisfies ( jTaruya, Nishimichi, Saito| 
land Hiramatsu|2009| ) 

^ dqn^{q;z)^C (65) 

where C = 0.18 in standard perturbation theory. SPT theory re- 
lies on a straightforward expansion of the set of cosmological hy- 
drodynamical equations and the approach has been repeatedly noted 
as being insufficiently accurate to model and describe the BAOs 
(peong and Komatsu 2006 , Taruya, Nishimichi, Saito and Hiramats^ 
[2009 , Nishimichi et al.2009 , Carlson, White and Padmanabhant2009r 
l^ruya, Nishimic hi and Sa ito 2010). In particular the amplitude of 
SPT predicts a monotonically increase with the wavenumber overes- 
timating the amplitude (Figure [2| with re spect to N-body simulations 
|Taruya, Nishimichi, Saito and Hiramatsu| ( |2009| . 

Note that we have divided the nonlinear power spectrum by a 
smooth linear power spectrum when constructing the ratio Rf [k) 
highlighting the scale dependence introduced by mode coupling. An 
alternative possibility would be to divide the nonlinear power spec- 
trum P^^ by a nonlinearly evolved smooth power spectrum Psmooth re- 
moving this scale dependence and providing a more detailed compar- 
ison of PT predictions against numerical simulations. We leave such a 
comparison to future work. 

Although we have not presented alternative approaches to PT 
there are a number of immediate extensions to this work that will be 
explored in upcoming papers. As future wide field surveys will cover 
both wide and deep regions of the sky we can use the SFB formal- 
ism as a tool to distinguish between different models for non-linear 
evolution of the matter density field. How do different theories affect 
the radial and tangential modes? How do different models affect the 
radialisation of information at various redshifts? 

We will compare different approaches to cosmological perturba- 
tion theory and investigate the power of the SFB formalism as a tool 
for model testing in forthcoming work. The role of other additional 
non-linearities, such as non-Gaussianity and modified galaxy cluster- 
ing bias models, will be presented in future work. 



4 RESULTS 

Following ' Rassat and Refreg"ier| ( |2012| | we construct the quantity 
Rf {k) to isolate the BAOs in the SFB formalism. The matter power 
spectrum, as defined in Eqn. ( [36| includes the physical effects of 
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Figure 2. Slice in 1-space showing Rf (k) for = 5 (top-panels) and ^ = 50 
(bottom-panels) in a wide and shallow survey of ro = 100/i~^Mpc (left- 
panels) as well as for a deep survey of ro = 1400/i~^Mpc (right-panels). The 
solid-lines denote the results of linear calculations and the dashed-lines show 
the total power-spectra computed by taking one-loop corrections into account. 
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baryons leading to the characteristic oscillations as seen in Fourier 
space (| Sunyaev and Zeldovich|1970[ [Peebles and Yu|1970| |Seo and| 
|Eisenstein,2 003 2007 ). 

In Figure (ij we construct slices of constant / through Rf {k) 
to investigate how non-linearities manifest themselves in the oscilla- 
tions. Rassat and Refregier ( 2012 ) used such slice plots to investigate 
the radialisation of information when varying levels of tangential and 
radial information is included in a survey. The radialisation of infor- 
mation can be investigated by notion that in the limit ro ^ oo we 
find 



lim Rf{k) - R^{k) = ^ 



\k) 



(66) 



Radialisation means that Rf {k) tends towards R^{k) in both phase 
and amplitude. This occurs as the tangential modes are attenuated 
due to mode-cancelling along the line of sight ( Rassat and Refregier 
|2012| ). The radialisation can be seen in Figures |[T) and ||2] (linear and 
un-redshifted lines) where the amplitude of the first peak is higher in 
the wide and deep survey than in the wide and shallow survey. Ad- 
ditionally the BAOs appear to only have a radial (k) dependence in 
surveys with a large radial parameter ro, as can be seen by the invari- 
ance the BAOs under a varying multipole /. The addition of redshift 
space distortions does not change this trend drastically though we do 
see more prominent radial and tangential dependence in Figure |[T) 
with the rate at which the BAOs radialise being naturally impacted 
due to mode-mixing, attenuation and peak shifts. The results appear 
to be in agreement with previous studies with percent level shifts in 
the peaks to smaller k and damping of the amplitude fNishimichi et al 
2007 , Nomura, Yamamoto and Nishimichi 2008 ; Smith, Scoccimarro 
and Sheth 2008 , Nomura, Yamamot o, Huetsi and Nishimichi||2009l 
Taruya, Nishimichi and Saito 20 10 ). As can be seen in Figure |[T) the 



BAOs seem to radialise, even in the presence of RSD, at large val- 
ues of the radius parameter ro, notably so for large multipoles / as we 
would anticipate (Rass at and Refregier|2012| ). 

The BAOs, as realised in SFB space, radialise as the survey size, 
ro, is increased. This corresponds to the amplitude and the phase of 
the BAOs tending towards the values as measured in the Fourier space 
ratio R^{k). As noted in Rassat and Refregier (2012), for a wide- 
field shallow survey the BAO will have smaller amplitudes and are 
spread across the (/, k) space. However, in a wide-field deep survey 
the unredshifted BAOs have larger amplitudes and are confined to the 
radial modes (e.g. we have lost the tangential, /-mode dependence). 
It was shown in ^ Rassat and Refregier] p012| > that the BAOs appear 
to radialise before the full SFB spectrum is able to and notably so at 
large / (Figures |3|4| ), one of the key motivations for implementing 
the SFB formalism. 

As can be seen, when the RSD are included then the radialisa- 
tion of information is damped and there is an increased dependence 
on the information contained in the tangential modes for genuine ra- 
dialisation to occur. The non-linear corrections, in contrast, radialise 
as expected but there is an increased dependence on tangential modes 
in the shallow surveys with the amplitude of the BAOs being amplified 
when higher /-modes are taken into account. 



5 CONCLUSION 

The baryon acoustic oscillations give rise to a characteristic signa- 
ture in the observed matter power spectrum that acts as a standard 
ruler. Unfortunately, the observed matter power spectrum is contam- 
inated and complicated by non-linearities arising from the evolution 
of density perturbations, galaxy clustering bias, redshift space distor- 
tions and survey specific systematic errors. In addition, upcoming fu- 
ture surveys will cover both large and deep areas of the sky demanding 
a formalism that simultaneously treats the both the spherical sky ge- 
ometry and the extended radial coverage. A spherical Fourier-Bessel 
basis was proposed as being a natural basis for random fields in this 
geometry. The recent study by [Rassat and Refregier|p012] ) was an ini- 
tial step into investigating the role of the SFB formalism in the study 
and analysis of the BAO. This study, however, did not go as far as 
including higher order contributions to the power spectrum that may 
impact the radialisation of information by introducing, for example, 
mode-mode couplings. The stability of this radialisation of informa- 
tion and the information content of tangential (/) and radial (k) modes 
for higher order physics is the key topic of interest. 

In this paper we have presented a short treatment of the effects 
of redshift space distortions and non-linear corrections to measure- 
ments of baryon acoustic oscillations in the spherical Fourier-Bessel 
expansion. In order to guide this investigation we have extended the 
formalism and techniques outlined in jRassat and Refregier| ( |2012|. In 
particular we have been able to use the procedure outlined in| Heav-| 
|ens and Taylor] ( ]1995] ) to construct a series expansion solution to 
model redshift space distortions. This solution was used to numeri- 
cally and analytically investigate the modulation to the angular SFB 
power spectrum. The qualitative behaviour of these corrections was 
outlined for surveys with varying levels of radial (/c-modes) and tan- 
gential (/-modes) information. It was seen that the RSD impact the 
radialisation of information through mode-mixing and the distinct sig- 
nature was investigated over a range of survey configurations. The ap- 
propriate formalism for partial- sky coverage was also introduced and 
outlined for which results will be presented in forthcoming work. 

The second main non-linearity investigated was one-loop correc- 
tions to the matter power spectrum arising from standard perturbation 
theory. A brief outline of perturbation theory methods was given and 
the basic equations for SPT introduced. The non-linear corrections 
and impact on radialisation of BAOs in the SFB was numerically in- 
vestigated. Possible parallel studies may include the analysis of dif- 
ferent perturbation theory methods and their ability to influence the 
radialisation of information as well as the possibility of distinguishing 
between theories for the non-linear corrections to the matter power 
spectrum and, consequentially, the BAO signal. Additionally we have 
neglected a detailed analysis of non-linear galaxy biasing models 

In this paper we have neglected other contributions to the power 
spectrum such as General Relativistic corrections, lensing terms, the 
role of non-linearities through more detailed studies, more complex 
treatments of galaxy biasing and more detailed modelling of the 
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Wide and Shallow Survey 




Figure 3. Ratio Rf {k) of SFB spectrum with and without the physical effects of 
baryons in (l,k) phase space for a wide and shallow survey of ro = 100/i~^Mpc 
using a Gaussian selection function. The baryonic wiggles are seen in both the 
radial (k) and tangential {t) directions. 



Wide and Deep Survey 




Figure 4. Ratio Rf {k) of SFB spectrum with and without the physical effects of 
baryons in (^,k) phase space for a wide and deep survey of ro = 1400/i~-'^Mpc 
using a Gaussian selection function. The baryonic wiggles are seen in both the 
radial (k) and tangential {t) directions. 



Wide and Shallow Survey 




Figure 5. Ratio Rf [k) of SFB spectrum with and without the physical effects of 
baryons in (l,k) phase space for a wide and shallow survey of ro = 100/i~-'^Mpc 
using a Gaussian selection function but with the inclusion of non-linear features as 
calculated in Standard Perturbation Theory. 



Wide and Deep Survey 




Figure 6. Ratio Rf [k) of SFB spectrum with and without the physical effects of 
baryons in (^,k) phase space for a wide and deep survey of ro = 1400/i~-'^Mpc 
using a Gaussian selection function but with the inclusion of non-linear features as 
calculated in Standard Perturbation Theory. 
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hydrodynamical and radiative processes involved in these processes 
{ [Guillet et al|2010[|Juszkiewicz, Hell wing and van de Weygaert|2012| . 
In addition we have not considered the role of systematic errors in a 
particular survey and how they may affect the radial and tangential 
measurements in different ways and consequentially complicate the 
analysis of BAOs. We leave this to a future paper but note that the 
results in |Asorey et al| ( |20l"2| will be of interest. 
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APPENDIX A: SPHERICAL BESSEL FUNCTION 



APPENDIX B: SPHERICAL HARMONICS 

The spherical-harmonics are complete and orthogonal on the surface 
of the sphere: 



I 



d(lYi^{d)Yi.^,{Q.') = 5fi,5, 



mm' 



(Bl) 



(B2) 



The overlap integrals of three spherical harmonics are given by the 
Gaunt integral which are expressed in terms of 3j symbols (denoted 
by matrices below): 
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APPENDIX C: 3J SYMBOLS 

The following orthogonality properties of 3j symbols were used to 
simplify various expressions. 
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POO ^ 

J r'^drji{kr)ji{k'r) = ^^^(^ - (^1) 

The first derivative of the spherical Bessel function can be expressed 
using the following recursion relation: 

'i{r) = [lji-i{r) - (/ + l)j^+i(r)] . (A2) 



Ji 



The second- and higher-order derivatives are deduced by successive 
application of the above expression: 

^^)=i(2/ + 3)(2/ + l)^^^^)- 

^JlMr) - ,07 , 10/ , oN -^^+2(0 • (A3) 



(2/ - 1)(2/ + 1)-^ 



(2/ + l)(2/ + 3)'' 



These expressions can be used to simply the kernels Ii'^\k, k') de- 
fined in Eqn.([33]) to express mode-mixing due to redshift- space dis- 
tortion. 



APPENDIX D: FINITE SURVEYS AND DISCRETE 
SPHERICAL BESSEL-FOURIER TRANSFORMATION AND 
PSEUDO-ClS 

Dl 3D Scalar fields 

Different types of boundary conditions are employed in the literature 
for finite surveys (IBinney and Quinn fT99T||Fisher et al|1995[[Heavens| 
land Taylor 1995 ). 

A natural choice for the boundary condition is to assume that the 
field vanishes at the boundary of the survey r = R leading to follow- 
ing condition on the radial modes that is determined by the zeros of 
the spherical Bessel functions ji (r): 



jiiqin) ji{kinR) 0; qin hnR. 



(Dl) 



The closure relation for spherical harmonics will take the follow- 
ing form: 
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^0 



dz z'^ji{kniz)ji{kniz) = ^[ji+i{qin)f S^n' ■ (D2) 



Which, in terms of the radial wavenumber, can be expressed as 
follows: 



^0 



dr r'^kinki>n'ji{hnr)ji'{kifn'r) 



(D3) 

The discrete spectrum is determined by the zeros of the spherical 
Bessel function. The normalisation coefficients are given by: 



1 W 



[klnjl + l{klnR)f ■ 



(D4) 

(D5) 
(D6) 



Tin 2 

The inverse and forward discrete SFB transforms are as follows 

^im(fan) = Tin J d^T ^ {r)jl{kr)Ylm{Cl)] 

^(r) = ^r^n^^m(A:)i^(/cr)y^^(Q). 

Imn 

The following expression is useful: 

r - 

^im(fon) = ^27t)^/^ J ^^k'^(^lr^,^k)Ylm{^k) (D7) 

In case of finite survey the 3D power- spectrum samples only discrete 
radial wave-numbers kin which is defined by the survey radius R: 



{'^lm{kln)^*l>^f{kifn') = P^i>^i>{kln)5w5mm'^n 



(D8) 



In addition to finite survey size, surveys often have a mask s{Q). The 
SFB transform of a masked field defines the convolved or Pseudo har- 
monics ^im(fon): 

^lm{kln) - \ -Tin f t'^ dv I dCl 

x[cl){r)s{n)]^{r)ji{kinr)Yirn{n)dn (D9) 

The convolved or P^^w Jo-harmonics are expressed in terms of all- sky 
harmonics ^im{kin) by the following expression: 



^l7n{kln) 



^ rVn'W{kln,kifr,')^lm{kvn' 



n' I'm' I" m" 

XSl''m"Ill'l'' 



(DIO) 



The kernel W{kin, ki^n') depends on selection function (/)(r): 

nR 

W{kin,ki^n')= / dr (t^{r)ji{kinr)ji{k'i.^.r) (Dll) 

The Pseudo-Ci s (PCLs) constructed from the convolved harmonics are 
a function of power spectrum of the angular mask Cf,, , normalisation 



coefficients Tin and the selection function ( 

Cl{kln) = {^lm{kln)^i^{kln)) 



Y^Y^v^ 2 j-ii'i" 

l^l^l^^i'^' 2/ + 1 



^ii'i" ( ^ ^' 







y.W^{klnM'r.')Cv{ki,n')Cin (D12) 

Notice that the PCls Ci{kin) are linear superposition of the power 
spectrum of underlying field Ci{kin). The mixing matrix Min^Vn' is 
given by: 



Cl(kln) - y^^Min^Vn'Cl{kvn' 



(D13) 



V n' 

where the mixing matrix is given by the following expression: 

If,,' ( S )V^(fc,„,W)Cf„. 

(D14) 

An unbiased estimates of the 3D power spectra can be written as: 



(D15) 



This is an extension of well known results for the projected surveys 
( [Hivon, Grski, Netterfield, Brill, Prunet and Hansen||2002j ). For low 
sky-coverage and small survey volumes the matrix Min^Vn' is ex- 
pected to be singular and binning of modes may be required. 

A different choice of boundary condition is often employed 
( [Fisher etalfT995l ): 



ji-i{k'inR) = 0; 
The normalisation constants in this case are given by: 



1 R^ 



[klnjl{klnR)f- 



(D16) 



(D17) 



The expressions for the mixing matrix derived above can still be used 
simply replacing the normalisation coefficients Tin- 

For discrete fields such as the galaxy distribution we can use 
the PCL approach if we replace the continuous function ^(r) 
with a sum of delta functions that peak at galaxy positions r^: 
^(j^) = Z^f^i^^^l^^ - rs); here N is the number of galax- 
ies. The SFB for such a discrete field is given by "^imik) = 
Z^fLi nmji{rskin)Yim{Cls)- Where the radial and angular position 
of galaxies are denoted by Ts = {rs,^s) = (^s, Os^(t)s) 



APPENDIX E: STANDARD PERTURBATION THEORY 

In the formalism outlined in section |3] any statistical observable can 
be computed to arbitrary order. Typically we are only interested in 
the second order corrections to the matter power spectrum though ex- 
pressions for higher-order corrections have been derived. One of the 
key issues regarding the inclusion is the computational costs required 
for these higher order corrections in part due to the high dimension- 
ality of the integrals, even after symmetry arguments have been taken 
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into account. The analytic expressions for the first corrections can be 
analytically derived [Makino, Sasaki and Suto| ( |1992] | 



1 A-^ P"^ 



42x^ + ^(x'-l)' (7x^ + 2) log 



12 2 
— - 158 + lOOx^ 



(El) 



1 + x 



1-x 



1 r 

P22{k) = — ^ y dxPy^{kx) j dfiPunik^ 

(3x + 7/i- lOx/i^)^ 
(l+x2-2x/i)^ 



1 + x2 - 2x/i) 
(E2) 
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